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Abstract 

Optical single-short measurement of the cross-correlation function between periodic sequences 
is demonstrated. The sequences are encoded into the broadband ultrashort phase-shaped pulses 
which are mixed in the nonlinear medium with additional amplitude-shaped narrowband pulse. 
The spectrum of the resulted four wave mixing signal is measured to provide the cross-correlation 
function. The high contrast between the values of cross-correlation and auto-correlation (the latter 
includes also the information of the sequence period) has potential to be employed in the optical 
implementation of CDMA communication protocol. 
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Modern communication technology, including the mobile 3G wide-band and Global Posi- 
tioning System standards, implemented in the rf-regime, extensively employs communication 
based on the Code Division Multiple Access (CDMA) protocol jj. It is known to be prefer- 
able over the Time Division and Frequency Division multiplexing strategies, where different 
time or frequency slots, respectively, are assigned to different users sharing the common 
communication channel. In CDMA the transmitter and receiver are assigned a unique code 
to recognize each other from the overall information flow. The codes are converted into the 
individual characteristics of different frequency components of the broadband spectrum and 
translated to the time-dependence of the waveform. This highly efficient spread spectrum 
multiplexing strategy is potentially capable to support a high number of users, given a large 
codes families with low correlation properties are available [2]. These code families, thus, 
must support sufficient number of users, while minimizing the interference between them. 
The state of the art code families incorporate periodic cyclically distinct pseudo- 

random sequences over the arbitrary alphabet. The codes have low correlating properties 
which imply minimal interference between them. This generalization to the extended alpha- 
bet dramatically improved the capacity and performance of the CDMA networks, while in 
the past only the binary codes were considered. 

It was recognized several decades ago that the developments of the communication tech- 
nology in the rf region are also applicable for the optical communication 6|, |7| . In the optical 
implementation the codes are spread over the broad bandwidth of the light pulses, which are 
phase-shaped in the frequency-domain |8|. Recently, the demand for the fast and reliable 
multiuser communication stimulated an extensive research activity in this field, focusing on 
different aspects of optical networking 

a a. 

Availability of a powerful technique to dis- 
criminate between different codes is critical to fully benefit from the advantages of CDMA 
0, [h]] . Here we demonstrate efficient all-optical experimental recognition of matching codes 
from the large family using single-shot measurement of their entire cross-correlation func- 
tion in the frequency- domain. The interference between the photo-induced optical pathways 
excited with two encoded broadband pulses and one amplitude-shaped narrowband pulse 
within the non-linear medium is employed to evaluate the cross-correlation function between 
the codes. We utilize the BOXCARS setup [ljj to obtain spatially-resolved background-free 
four-wave-mixing spectrum which is formally equivalent to the cross-correlation function be- 
tween the codes of interest. Our approach enables the high-contrast discrimination between 
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the periodic low-correlating codes based on identification of the constructive interference 
corresponding to the pulses shaped with similar codes and destructive interference corre- 
sponding to t he p ulses shaped with different codes. As opposed to the previously developed 
techniques (9), [ly], the presented setup enables single-shot measurement of the entire cross- 
correlation function, which in addition provides usefull information on the unknown period 
of the code without additional computational expences. 

We consider the families of periodic cyclically distinct pseudo-random sequences Sk(t) 
over the alphabet of integers of size p [0 < Sfc(t) < p — 1] with arithmetics modulo p, where 
p >2. The cross-correlation function between two sequences is defined by 

c(k,r,r) = J2^ t+T) - sAt \ (1) 

i=0 

where k and r are indexes of the sequences in the set of size M, < t < L — 1, L - is the 
period of the sequence, and 

Up = e i2n / p - is pth root of unity. The commonly used merit 
factor of the codes performance is the maximal value of the correlation C max =max{c(k, r, r)} 
excluding the cases where simultaneously k = r and r=0, for which, trivially, C max = L. 
We thus demand C max <C L to ensure the non-intercepted communication with the absence 



of interference. The nested chain 



ow-correlation codes families S(0) C S'(l) C S(2) C . . . 



obey the desired requirements 



5j. For a given period L, the code family of higher 
order incorporate more code sequences, while the corresponding value of C max increases. 
However, as the code period L increases, the number of available codes increases, while the 
value of C max decreases. The C max values of the S'-families have a well defined bounds and 
distribution, and thus can be appropriately tailored to the implementation needs by the 
choice of L and the order of the family. 

In the optical domain, the common implementation of the CDMA network communication 
involves broadband ultra-fast shaped laser pulses . The high popularity of this benchmark 
system is mainly due to the Weiner- Heritage design of the pulse shaping apparatus that is 
used to assign the codes to the laser pulses carrying the transmitted data. The 4f optical 
setup of the pulse shaper incorporates pixelated liquid crystal Spatial Light Modulator 
(SLM) in its focal plane, where each spectral component of the laser electric field E(u) = 
\E(u) \ e %q> ^ is addressed individually. Here, - is the spectral amplitude at frequency 

u, and - is the corresponding phase. The phase-shaping is implemented via control 

over the voltage applied to different pixels of liquid crystal, which induces different phase 



retardation of the transmitted light. We assume that spectral segment of Auj is spread along 
each pixel. The similarity between the representation of integer code sequences s k (r) in the 
complex plane and the phase factor of the electric field is used to convert the codes into the 
relative phases of the broadband pulse spectral components: $(c<j) = (27r/p)s k (u>). Once the 
encoding process is completed at the transmitter site, measurement of the cross-correlation 
function c(k, r, r) is the important step that should be performed at the receiver or router 
site. 

Naturally, the single-shot all-optical cross-correlation measurement is the ideal way to 
effectively discriminate between different codes and identify their period. Translation of 
the codes into the phases of the pulse spectral components allows one to take advantage of 
the interference between multiple photo-induced multi-photon pathways excited in the non- 
linear medium by the broadband pulse. The cross-correlation function c(k, r, r) between 
the two codes s k and s r converted into the phases $ k (u) and $ r (c<j) of the pulses with the 
electric fields E k (oo) and E r (u) obtains the form of 

/OO /»00 
E k (u)E* r (u; + A)dw = / \E k (u)\ ■ \E r (u + A)\S^~^ u+A ^du. (2) 
-OO J — OO 

Neglecting the spectral envelope profiles of \E kjr {uS)\, r, A) is formally equivalent to 
c(k,r,r), differing from it by A being continuous variable representing detuning between 
spectral components of E k {uj) and E r {uj + A), while r is a discrete integer variable repre- 
senting the relative delay between the two codes. 

Physically, r, A) interferes all the possible two-photon pathways of one absorbed 
photon and one emitted photon which spectral components are detuned by A within the 
Raman-type process. Suppose that variable u is divided to small bins of size Auj and each 
bin is associated with specific value of t - the index of the code component, so that the 
phase $(cj) of the electric field E(uj) is constant across each bin. For simplicity, we analize 
first the case where the spectral envelope is extremely slow and can be approximated to 
be constant: \E k ^ r {uj)\ ~ 1, Similar to c(k,r,r), r, A) obtains its maximal values when 
A = nAuoL if k = r, as the interference between all the corresponding two-photon pathways 
is fully constructive. Here n = 0, 1, 2 . . . and L is the period of the code. For A ^ nAuoL 
or k 7^ r the interference between the two-photon pathways is destructive and the value of 
r, A) <^ k, nAujL). Overall, the profile of r, A) has peaks corresponding to the 
cases of the constructive interference that appear with the spacing associated with the period 
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of AuL. If the interference is destructive, for the case of flat spectral profiles of E k)r {u), the 
bounds on the value of £(k, r, A) are obtained from the results of analysis of c(k, r, r) within 

nnn 

the framework of Number Theory [3J, |4J, |5| . For the proper choice of the codes family, these 
interference is kept sufficiently low and enables clear discrimination between the matching 
[k = r] and non-matching [k ^ r] codes. 

Experimentally, multiple two-photon pathways are interfered within the non-linear 
medium. However, due to its Raman nature, the presented two-photon photo-excitation is 
not applicable for the single-shot detection. The latter is possible employing the four-wave 
mixing process within the medium having third-order susceptibility x ■ The interaction 
includes absorption of two-photons and emission of one photon provided by the excitation 
electric fields of the lasers, while the fourth emitted photon generated by the non-linear 
polarization of the medium is detected: 

/oo poo 
E pump (Lu')duj' E k {u")E*{u" + u'-u)du". (3) 
-oo J — oo 

To avoid the broadening of the interference pattern created by the pulses Ek{uj) and E r (u) 
due to the presence of the additional broadband pump laser pulse E pump (uj), the latter must 
be of narrow bandwidth. Let us assume that 



pump 



E pump (uj')5{ujQ - u)')du'. (4) 



Then we obtain 



/oo 
E k (u)E*(u + A)du (5) 
-oo 

/oo 
\E k (u)\ \E* r (uj + A) I e ^H-WA)]^ 
-oo 

Therefore, we establish the recognition of matching codes and measuring their period upon 
the measurement of the codes cross-correlation function in the frequency domain, which is 
equivalent to the measurement of the spectrum of E sig {uo). In the case of the finite bandwidth 
of E pump the measured spectrum of E sig {ui) represents the cross-correlation function £ (k, r, A) 
convoluted with spectral profile of E pump . 

Our experimental setup utilizes ultrashort pulses of 65fs (fwhm=15nm) produced by 
Ti: Sapphire oscillator at 800 nm that were regeneratively amplified at lKHz rate and splited 
into three beams. Two beams were passed through the phase-shaping 4f optical set-up incor- 
porated 640-pixel liquid- crystal Spatial Light Modulator which was used to assign different 



codes for each beam independently. The corresponding experimental shaping resolution 
was Au;=0.11 and 0.136 nm/pixel respectively. The third beam was amplitude-shaped us- 
ing similar 4f shaping set-up to yield narrow [fwhm=0.9nm] spectrum. The resulted three 
beams were focused at zero delay within the BOXCARS geometrical configuration onto the 
lmm-wide piece of fused silica. The resulted spatially isolated signal was collected and its 
spectrum was measured by CCD camera coupled to the spectrometer with resolution of 
O.lnm. It worth noting that alternatively to the presented BOXCARS geometry, polar- 
ization gating or self-diffraction geometries can be implemented within the same general 
formalism. 

We have used the quaternary [p=4] codes which belong to the nested families 5(0) and 
5(1) and generated using corresponding linear shift registers constructed upon the prim- 
itive polynomial. The chosen period of the codes was L = 2 l — 1=7, where Z=3 is the 
order of the generating polynomial. The size of the family 5(0) is M = L+2=9 with 
\C max \=^/LTl+l~3.8. The size of the family 5(1) is M = (L + 2)(L + 1)=64. Formally, 
the corresponding value of C max for this case \C max \—2y/L + 1+1—6.66, however, for the 
specific case of 1=3, occurrence of this value is zero and the effective \C max \=5. 

On Figure [1] the typical spectrum of E sig (uj) is shown for the cases of k = r, and k ^ 
r. The very pronounced peaks observed in the spectrum correspond to the constructive 
interference between the photoinduced pathways associated with the identical codes, where 
k = r. The peaks are separated by A^ = nAuL, where n - is the integer parameter. This is 
in high contrast to the spectrum obtained for the case of the destructive interference among 
the pathways for the case of different codes, i.e. k ^ r. For the peaks exceeding the certain 
threshold, the distances between the peaks is used for the determination of the code period 
corresponding to different values of n. The corresponding results for the codes belong to 
the families 5(0) and 5(1) are shown on Figure [2] in the form of 2D histogram, where x- 
axis corresponds to the programmed code period and the y-axis - to the retrieved period. 
We have successively recognized periods of 98% out of total 432 different periodic codes. 
As it is seen from Figure [2] the origin of the 2% of the false recognized codes comes from 
the shortest and longest periods examined. For the shortest period, the factor limiting the 
recognition power is the actual width of the probe pulse, while its convolution smears the 
spectral features and broadens the spectrum. On the side of the longest periods, the total 
spectral width of pulses E k ^ r (uj) becomes important. If the pulse bandwidth is not broad 
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enough, there is not enough constructively interfering pathways provided and the contrast 
between the cases oik — r and k ^ r is not sufficient for the proper recognition. 

The codes with different periods are easily distinguished form the measurement of E sig {uj) 
due to the lack of the constructive interference to the amplitude of the emitted signal. For 
the experimental discrimination between different codes with the same period it is enough 
to measure the signal at specific frequencies E sig (ujo + mAj,), where m = 0,1,.... We 
present the corresponding results for the code families S(0) and 5*(1) on the 2D color-map 
plot of Figure [3] (a) and (b) respectively. Here the color represents the normalized value 
°f YlmiEsigiyjQ + mAi) for the different pairs of codes k [x-axis] and r [y-axis], so that 
the auto-correlation cases appear on diagonal. As the signal corresponding to r = k is 
sufficiently higher than for r 7^ k, we obtain a very efficient recognition of the 100% of codes 
belong to S(0) with the signal threshold of 90% [see panel (a)]. For the codes of the ^(l) 
family the value of C max is higher then that of S(0) and, hence, the threshold for the code 
recognition must be also higher. We obtain 0.6% of experimental false identification using 
threshold of 95% in this case [see panel (b)] out of the whole set of the cross-correlation 
measurements. Our experimental results were verified with the numerical evaluation of 
the spectrum of E sig (uj) considering experimental parameters as pulse-shaping resolution, 
bandwidth of the pump pulse, etc. The same fraction of the false identified signals was 
obtained. We have also performed numerical simulation of the signal considering the pulses 
with broader bandwidth [fwhm=25nm, pulse duration 45fs] , our experimental pulse-shaping 
resolution, and bandwidth of the pump pulse of fwhm=0.8nm for the S(l) family with 1=4. 
Here L=15, M=255, while \C max \=9. For this case we obtain that 100% of codes can be 
recognized applying the 90% threshold. 

In conclusion, we have experimentally demonstrated the all-optical single-shot frequency- 
domain cross-correlation measurements for the efficient matching and simultaneously mea- 
suring the period of low-correlating codes. The presented technique has potential to be 
implemented in the optical phase-shaped spread-spectrum code division multiple access com- 
munication network. 
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FIG. 1: Typical experimental spectrum of E s i g (u) for the cases of identical [k = r, auto-correlation, 
black line] and different codes [k ^ r, cross-correlation, red line]. The distance between the peaks 
for the case of k = r corresponds to the period Ax,. 




Programmed period, (Acq/.) 



FIG. 2: 2D histogram of measured vs programmed code period Ax for 432 different codes belong to 
family S(0) and S(l). The distribution along the right axis reflects the signal-to-noise experimental 
ratio, which results in the 98% true matching of the code period. 
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Code index [family S(0)], k Code index [family S(1 )], k 

FIG. 3: Color-map plot representing the maximal cross-correlation values E sig (ujQ + mA^), 
m=0,l,. . . for different codes, (a) - results for the family S(0), 1=3. 100% of the codes was 
efficiently matched based upon the cross-correlation measurements applying the threshold of 10%. 
(b) - family S(l), 1=3. Applying the threshold of 95% we obtain that 0.6% out of whole set of 
measurements led to the false identification. This result agrees with our numerical evaluation 
considering the experimental parameters. 
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